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ABSTRACT 

A recent study of ultra-compact dwarf galaxies (UCDs) in the Virgo cluster revealed that some 
of them show faint envelopes and have measured mass-to-light ratios of 5 and larger, which 
can not be explained by simple population synthesis models. It is believed that this proves that 
some of the UCDs must possess a dark matter halo and may therefore be stripped nuclei of 
dwarf ellipticals rather than merged star cluster complexes. 

Using an efficient N-body method we investigate if a close passage of a UCD through 
the central region of the host galaxy is able to enhance the measured mass-to-light ratio by 
tidal forces leaving the satellite slightly out of virial equilibrium and thereby leading to an 
overestimation of its virial mass. 

We find this to be possible and discuss the general problem of measuring dynamical 
masses for objects that are probably interacting with their hosts. 

Key words: galaxies: dwarfs - galaxies: interactions - galaxies: kinematics and dynamics - 
methods: N-body simulations 



1 INTRODUCTION 

Hase gan et alJ 0005) investigated ultra-compact dwarf galaxies 
(UCDs) around M87, the central galaxy of the Virgo cluster. By 
measuring the surface brightness profiles and assessing the pro- 
jected velocity dispersion of the objects they concluded that some 
UCDs of their sample have mass-to-light ratios of the order 5- 
9. Furthermore, they find that some of the UCDs show faint en- 
velopes. This supports the notion that these UCDs may be stripped 
nuclei of dwarf ellipticals. 

UCDs were first disc overed by Hilke^ ll998l) : 
Hilke rThifante & Richtlerl J 19991) during a study of globular 
clusters and dwarf galaxies around the central galaxy in the Fornax 
cluster. These objects are compact with effective radii of about 
15-25 pc and a hundred to several hundred pc in extension, and 
they are massive with masses of a few 10 6 up to a few 10 7 Mq. 

There are several theories about the origin of UCDs: (I) They 
could be the most luminous end of the distribution function of very 
massive globular clusters (Hilk er et alJI 199 9 : Mieske et al. 2002; 
iDirsch et alJ 120031). (II) They could be the remnant s of st ripped 
dwarf ellipticals jBekki. Couch & Drinkwate3l200lt iBekki et all 
l2003t iDe Propris et alJl2005l) . In this 'threshing' scenario a nu- 
cleated dwarf elliptical looses its envelope and most of its dark 
matter content due to tidal interaction with the host galaxy, such 
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that only the 'naked' nucleus remains. (Ill) They could be amalga- 
mated young massive star clusters formed in a star cluster complex 
during the star-burst caused by the tidal perturbation and possible 
disruption of a gas-rich galaxy iFellhauer & KroupjEo02ll2005h . 
This scenario is well-established theoretically and was first pro- 
posed by Kroupa <1998h . It must have been profusively active dur- 
ing the early hierarchical structure formation epoch when gas-rich 
substructures merged to the present-day major galaxies. 

It is still under debate whether the UCDs, which fill the gap 
between globular clusters and dwarf galaxies, follow the fundamen- 
tal plane (effective radius / velocity dispersion - total luminosity) 
relation for globular clusters or for dwarf gala xies lHasegan et alJ 
l2005UEvstigneeva Gregg & DrinkwateJ2005h . While, in the M v - 
gn-plane (see e.g. lHasegan^HrtT2o"o5l their Fig. 7), they lie closer 
to the globular clusters, their relation seems to rather follow the one 
of dwarf galaxies. Furthermore they occupy exactly the space be- 
tween dwarf ellipticals and their nuclei. This seems to point to for- 
mation theory (II). But the surface brightness profiles of the bright 
UCDs in Fornax are much more extended (i. e. larger effective radii) 
compared with nuclei of dEs foe Propris et al. 2005). 

On the other hand, M araston et al] 120 04) found an interme- 
diate age object (W3, age 300-500 Myr) in the merger remnant 
galaxy NGC 7252. The mass (M = 8 • 10 7 M Q ), size (r cff = 
17.5 pc, and velocity dispersion of 45 kms -1 ) strongly suggest 
this object to be a UCD rather than a globular cluster. The age of 
this object, which corresponds to the time elapsed since the ma- 
jor interaction, unambiguously shows that this object can not be a 
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Table 1. Table of our initial model parameters. The columns denote the initial mass of our 
model (M p i = M; n ;), the scale length (Plummer radius, i? p i; analytical & measured as 
described in Sect.|3}, the characteristic crossing time (T Gr ), the central projected (line-of-sight) 
velocity dispersion (cto, p ; analytical & measured as described in Sect.[3j, and the scaling factor 
(A) to compute the virial mass (see Eq. I1 II . 



Mass [M ] 


Rpi [pc] 


measured 


T cr [Myr] 


co,p [kms x ] 


measured 


.4 


10 7 


25 


23.7 


3.69 


15.92 


16.20 


1608 


10 7 


50 


46.6 


10.43 


11.25 


11.71 


1565 


10 7 


100 


94.6 


29.50 


7.96 


8.43 


1487 


10 7 


250 


221.0 


116.61 


5.03 


6.30 


1140 


10 s 


25 


25.2 


1.17 


50.33 


50.2 


1573 


10 x 


50 


47.4 


3.30 


35.59 


36.8 


1558 


10 x 


100 


92.8 


9.33 


25.16 


26.9 


1489 


10 s 


250 


213.0 


36.88 


15.92 


19.91 


1184 



stripped nucleus of a dwarf elliptical. A dwarf elliptic al can not be 
strippe d in a time interval of only 500 Myr. iFellhauer & Kroup3 
(2005) showed that W3 could be the merger object of a massive 
star cluster complex which was formed during the interaction. The 
subsequent merging of star clusters forms an object with properties 
similar to those of W3. The evolution of the simulated super cluster 
shows t hat it trans forms into a UC D such as those fo und in Fornax 
( Hilke r et alll99Sl : IPhillipps et all200lh Abell 1689 <Mieske et alJ 
l2004l) , an d Virgo Masegan et alJl2005t lEvstigneeva et alJl2005h . 
Moreover jFeHhaue^^Kroup^T2005l) showed that, due to its high 
mass, the object was able to retain an envelope of bound stars 
which initially were expelled from the individual clusters during 
the merger process. Thus, an envelope around a UCD is not a proof 
of its cosmological origin as a dwarf elliptical. Still the puzzle of 
high mass-to-light ratios of some of the UCDs in Virgo remains. 
These mass-to-light ratios (5-9) are found for the UCDs lying clos- 
est to the host galaxy. This suggests that deviations from virial equi- 
librium may play a role. 

The reasoning is that dwarf galaxies are on radial rather than 
circular orbits, if the UCDs formed as star-cluster complexes during 
the merging of major gas-rich substructures. This allows the satel- 
lites to pass close to the galactic centre. Tidal forces in the central 
region of the host galaxy are rather strong. Hence, the dwarf object 
looses stars and subsequently departs from virial equilibrium. Some 
of the lost stars (stars which are no longer bound to the object) do 
not immediately leave the object or its vicinity but disperse slowly 
along the orbit. Thus, a line-of-sight velocity dispersion measure- 
ment may be contaminated by these unbound stars which inflate 
the velocity dispersion. Furthermore, the gravitational shock of the 
central passage leads to an expansion of the dwarf galaxy which 
is later reversed again. But still this expansion through tidal heat- 
ing may result in a measurable increase of the core radius once the 
dwarf object is again outside the host galaxy. Both of these effects 
may lead to an overestimation of the measured dynamical mass re- 
sulting in a higher mass-to-light ratio. These effects are not new to 
the astronomi cal comm unity and are studied intensively by various 
autho rs (e.g . lKroupalll997l) for dwarf spheroidals, |Ma^er^yilJ 
2001, 2002) for dwarf discs and dwarf spheroidals or l Bek ki et alJ 
2003) for dwarf ellipticals). With this paper we want to extend these 
studies onto massive (compared to globular clusters) and compact 
(compared to other dwarf galaxies) objects like UCDs. 

It is definitely clear that these effects must be very strong if 
an object comes close to the galactic centre. But on the other hand 
these very close passages may be highly unlikely. By means of nu- 



merical simulations we intend to find out which sets of orbits allow 
an enhanced mass-to-light ratio, i.e. how close the UCD has to ap- 
proach to the centre of its host galaxy, for which we choose the pa- 
rameters of M87 to account for the UCDs with high mass-to-light 
ratios found around the central galaxy of Virgo. 



2 SETUP 

We model the parent galaxy as an analytical potential because on 
one single passage dynamical friction for a dwarf galaxy of mass 
M ?C 10 s Mq affects the orbit t o at most 2-3 per ce nt, estimated 
us ing Chandrasekhar's form ula I Chandrasek haJ 1 943h as described 
in lPortegies Zwart & McMillanl 120021) . We also look for only one 
central passage because unbound stars disperse along the orbit and 
are lost at subsequent passages. Some of them might still be around 
the object after the second or third central passage but definitely not 
for dozens of orbits. 

As a model for the analytical potential we choose the param- 
eters for M87, the central galaxy in the Virgo cluster, consisting 
of a NFW-profile for the dark halo, a Hernquist profile (H) for the 
visibl e matter (stars) and a central super-m assive black hole (BH) 
I McLaughtir] 1 1999t IVesperini et"ai]l2003t |Pi Matteo etai]|2003l) . 
The density profile of the host galaxy is 



Ptot(r) 



= pNPw(r) + Ph(v) 

PO.NPW f Si NFW 

1 + — ^~ 



+ 



Mh r a ,H 



27T r (r s , H + r) 



with the following parameters 

3.17- 1C 

560 kpc 



P0,NFW 
r s ,NFW 

M H 

J~s,H 

M BH 



3.17- lO" 4 M /pc 3 , 



5.1 kpc, 
3- 10 9 M . 



(1) 



(2) 
(3) 
(4) 
(5) 
(6) 



The above parameters denote from top to bottom the characteristic 
density and the scale-length of the NFW-profile, the total mass and 
the scale-length of the Hernquist-profile and finally the mass of the 

central super-massive black hole. 

We model the UCD as a Plummer- sphere <Plummerll91ll) in 
the numerical realisation described by lAarseth. Henon & W ielen 
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Table 2. List of minimum distances and the corresponding tangential ve- 
locities at the start of the simulation. 



D m j n [pc] titan [kms *] 





0.0 


50 


6.0 


100 


10.9 


150 


15.9 


250 


25.4 


500 


48.1 


1000 


89.3 


1500 


126.0 


2000 


158.9 



119741) . 



3M p i 



1 + 



R 2 



-5/2 



(7) 



with varying scale-leng ths, R p i, and initial masses, M p i, because 
iDe Propris et alj 120051) found a Plummer-profile to fit four out of 
five UCDs in the Fornax cluster quite well. Furthermore, the Plum- 
mer model is analytically simple and the Plummer radius is not 
only the scale-length of the model but also its half-light radius (the 
projected radius from within which half of the light of the object is 
emitted). 

The masses of UCDs range between several million Mq up 
to several tens of millions. So we choose for our models 10 7 and 
10 s M© as initial masses (M p i = M- ln {). We vary the initial scale- 
length of our models to be R p \ = 25, 50, 100 and 250 pc to deter- 
mine the influence of the concentration of the objects. The cut-off 
radius (i?n m ; the radius where we truncate the Plummer distribu- 
tion of our object) of all our models was kept constant at 500 pc 
which is larger than the initial tidal radius. A detailed list of our 
model parameters can be found in Table Q The Plummer spheres 
are modelled using 10 6 particles. The object is set-up and inte- 
grated in isolation until equilibrium, as measure d by the constancy 
of the 90 % Lagrangian radius, is reached iKroupal 19971) . 

The UCD model is then placed at a distance of 10 kpc to the 
centre of the host galaxy with no radial velocity. This means that the 
effect of the central passage is the most harmful possible because 
it is the slowest possible. The faster the passage would be (if the 
satellite was to start further out) the less tidal influence the central 
passage would have. No UCD is found closer than 10 kpc from the 
centre of its host and therefore it acts as a minimum apogalacticon, 
which has the strongest tidal effect possible. To vary the minimum 
distance (perigalacticon) we give our models different tangential 
velocities which are listed in Tab. [2] We assume the problem to be 
spherically symmetric, so we are able to place the trajectory of our 
model in the z-y-plane of our simulation area without restricting 
the problem. 

We use the particle-mesh code SUPERB OX to carry out the 
simulations. SUPERB OX has high-resolution sub-grids which stay 
focused on the core of the dwarf object while it is moving through 
the host galaxy. The resolution of the innermost grid contain- 
ing th e core is 3 pc. For a detailed description of the code see 
Fellha ueretail 120001. 



3 RESULTS 

We carried out a parameter survey of 76 simulations covering dif- 
ferent dwarf galaxy objects and different orbits. The parameter 
range is shown in TablesQand|2| For each parameter set one sim- 
ulation over two central passages is carried out. For the determina- 
tion of our results we look at the satellite when it reaches apogalac- 
ticon again after the first central passage. 



3.1 The analytical method 

The surface density profile of the satellite is fitted by a Plummer 
profile of the form 



E(r) = So 11+-^-) , (8) 

by applying a non-linear least-squares Marquardt-Levenberg algo- 
rithm. From this procedure we take the fitted Plummer radius for 
our determination of the virial mass. It can be shown that the Plum- 
mer radius exactly coincides with the half-light (projected half- 
mass) radius of a Plummer sphere. 

Furthermore we determine the line-of-sight velocity disper- 
sion profile. In cases where it is possible, i.e. the profile is not too 
contaminated by unbound stars, we again fit the Plummer profile 
for the line-of-sight velocity dispersion, 



= "Up 



1 + 



*8. 



-1/4 



(9) 



The projected velocity dispersion, as well as the surface den- 
sity, is measured along all three Cartesian coordinates in logarith- 
mically spaced, concentric rings centred on the object. For the mea- 
surement all stars within a certain distance i? max in front and be- 
hind the centre of the object are taken into account. The procedure 
is done along all three Cartesian axes and we take the arithmetic 
mean value (i.e. a mean profile), because we do not know which 
orientation our objects have with respect to the observer. Therefore 
effects may be very strong measuring along the trajectory of the 
object but almost not visible in the perpendicular direction. Any 
random direction is likely to yield an intermediate result. 

For the surface density we always fit a Plummer profile even if 
the object is completely destroyed and does not follow a Plummer 
distribution at all. But in most of our models the remaining object is 
still fairly well represented by a Plummer profile even if it is on the 
way to complete destruction. In cases where we are not able to fit 
the velocity dispersion profile with a Plummer profile (see e.g. first 
panel of Fig. [5} we take an average of all measured line-of-sight 
velocity dispersion values within the innermost 10 pc in projected 
radius to determine the central line-of-sight velocity dispersion. 
We now determ ine the virial mass taking the formula from 
2005 



Hasegan 
King, 7l96i 



M vir = 



et alJ 
[ST 



which is based on the theoretical work of 



9 V B 2 



(10) 



where R c is the core radius and <ro,p is the central value of the pro- 
jected velocity dispersion. The parameters v, a, and p are depen- 
dent on the kind of King model one wants to fit. Because we are 
using Plummer spheres we accumulate these parameters together 
with the constants into one single parameter A, which we determine 
for the isolated Plummer sphere before we start the simulation, 
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Figure 1. The ratio between virial and bound mass (M v - n - / M^,) for our ob- 
jects, measured at t = 60 Myr, the point where the object is again most 
distant to the centre of the host galaxy, against the closest distance from the 
centre of the host (D m i n ). First panel: Objects with initial mass of 10 7 Mq. 
Second panel: Objects with initial mass of 10 s Mq. Tri-pointed stars (on- 
line: black) denote objects with initial Plummer radii of R p i = 25 pc, 
crosses (online: red) have R p \ = 50 pc, five-pointed stars (online: green) 
have R p i = 100 pc, and six-pointed stars (online: blue) i? p i = 250 pc. 

where the values of the scaling factor A can be found in Table Q 
As one can see these values are lower for higher masses and in- 
crease for larger scale-lengths. Nevertheless we use the same A de- 
termined for the isolated model for all final models stemming from 
this initial model, even if the mass-loss is significant. This may lead 
to a slight underestimation of the final virial mass. 

On the other hand SUPERBOX calculates the number of bound 
particles (energy below zero) at each time-step. Therefore we know 
the real bound mass Mb. 

With these two masses we can determine a virial-mass-to- 
bound-mass ratio (M v i r /M^ or short M/M) which can then be mul- 
tiplied by a 'normal' mass-to-light ratio for a population of stars 
with the determined age and metallicity and without dark matter to 
obtain the dynamically measured M/L-ratio. 

This M/M-ratio is plotted in Fig.Q when the satellites have 
reached their apogalactica again. At this point the satellites are al- 



Figure 2. Parameter space of our simulations; Plotted are the Plummer ra- 
dius of our models, iZpi, against the closest distance to the centre of the host 
galaxy, D mln . Capital D (online: red) denotes the satellite gets completely 
dissolved during the first passage. Small d (online: magenta) denotes satel- 
lites which do not survive the next central passage. Small s denotes satellites 
which survive and form bound objects even after the next central passage. 
Capital E marks simulations where the M/M-ratio is enhanced if the central 
values are directly measured. Small e marks simulations which only show 
enhanced M/M-ratios if the observational method is applied (see main text 
for explanation). Surviving satellites with enhanced M/M-ratios are plot- 
ted green online. Finally, n marks simulations which show no enhancement 
measured with either method (online: blue). First panel is for objects with 
initial mass of 10 7 Mq, second panel for 10 8 Mq. 



most back to virial equilibrium again if they are not on the way to 
complete dissolution. 

The results show clearly that only objects which are not com- 
pact and/or have a low mass can be influenced enough by one cen- 
tral passage to enhance the M/M to account for the high mass-to- 
light ratios found in UCDs. But these satellites will also not survive 
the next central passage or have not survived the first one at all. 

The best representation of real UCDs is the model with an 
initial mass of 10 7 Mq and an initial Plummer radius of 25 pc. The 
results for this model are shown as tri-pointed stars in the first panel 
of Fig.Qand the lowest line in the first panel of Fig.|2| The results 
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Figure 3. Line-of-sight velocity dispersion profiles and the Gaussian fit values for three of our satellites. Three-, four- and five pointed stars denote the 
line-of-sight velocity dispersion values measured in concentric rings around the satellite along the x-, y and z-axis, respectively. Filled squares (online: red) 
are the mean values of all directions. Open symbols (online: magenta) on the right denote the central line-of-sight velocity dispersion values derived with the 
observational method (see Fig.|5|. First panel: This object gets completely dissolved and shows an enhancement of the velocity dispersion in either method. 
But this destroyed satellite can not account for the UCDs in Virgo (see Fig.|4j. Second panel: Object which shows an enhancement only if the central velocity 
dispersion is derived with the observational method. The open symbols denoting the observational values in the x- and y-direction are clearly above the 
'real' central value. Third: Simulation of a very massive satellite which shows no clear enhancement in both methods. 




Figure 4. Contour plots of the three models of Fig.|3]in the orbital plane (x-y-plane). The centre of the host galaxy is at 0-0 kpc. The orbit of the object 
is plotted as the thick (online: green) line. Contours have magnitudo spacing (assumed M/L = 1.0 Mq/Lq) and the outermost contour corresponds to 
28 mag arcsec (the two brightest contours in each panel are plotted red online). Almost filled black areas do not denote high surface brightness but areas 
with very faint contributions (at the given resolution the brightness fluctuates from pixel to pixel at the lowest level). First panel: This object is clearly 
dissolved; it shows no core anymore and is out of virial equilibrium. Second: Surviving object which shows faint tidal tails along the orbit which enhance 
the measured velocity dispersion. Third: This object shows only a few stars which are unbound and spread (black areas). The small fraction of these stars 
are not able to enhance the velocity dispersion. 



show clearly that either the satellite gets completely dissolved if 
D m in is closer than 50-100 pc or shows no deviation in M/M at all 
(Fig.0. 

In Fig. |2| we distinguish the models according to their sur- 
vival of the passage and if they show enhanced M/M-ratios. The 
panels show clearly that only dissolved or almost dissolved objects 
show enhanced M/M-ratios. All surviving objects show no devia- 
tion from virial equilibrium. 



3.2 The 'observational' method 

The results in Sect. 13. ll are based on the exact knowledge of the 
theoretical central line-of-sight velocity dispersion. Observers on 
the other hand do not have this information. The usual way to de- 



termine the velocity dispersion of a marginally resolved object is to 
place a slit on the object and obtain a spectrum. From this spectral 
information one chooses one or two spectral lines and fits a tem- 
plate spectrum convolved with the instrumental line- width function 
and a Gaussian velocity distribution. This Gaussian velocity distri- 
bution measures the line-of-sight velocity dispersion of the whole 
object. Based on assumed theoretical models this value is then cor- 
rected to the central line-of-sight velocity dispersion. In the case of 
the simple Plummer model one has to multiply the result for the 
whole object by a factor of 1.25. Thus, the central velocity disper- 



0,p 



3ttGM p i 
64R pl ' 



(12) 
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Figure 5. Same simulations as in Fig.|3] The plots show the velocity distributions along the x-axis. Histogram (online: black) are the data-points. Solid line 
(online: red) is the single Gaussian fit, dotted lines (online: green) are the three Gaussians to correct for fast and slow unbound particles in front and behind 
the object and the dashed line (online: magenta) is the sum of these three Gaussians, which is an excellent reproduction of the data and is barely visible 
by being masked by the data. The three Gaussians always fit the data best. While in the first panel the object is almost dissolved and the Gaussians of the 
unbound particles are stronger than the bound one (resulting in a M/M-ratio of 20 in x and 5 in j/-direction), the second panel shows a bound object with 
enhanced M/M-ratio due to the small contribution of the unbound particles (M/M-ratio is 2.5 and 2.2). In the last panel the contribution of the unbound 
particles is too low to enhance the M/M-ratio (M/M = 1.2 in both directions and within the uncertainties of the measurement). 



while the projected velocity dispersion integrated over the whole 
Plummer sphere is 



7 Pl 
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p 1 Jo 



— — / 2nr E(r )o>(r )dr , 



3ttGM p i 
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(13) 



= 1.25. 



(14) 



so that 

pi 
pi 

°obs,p 

<7q p denotes the central line-of-sight velocity dispersion and cr^ s 
the weighted line-of-sight velocity dispersion integrated over the 
whole object. 

As one can see in Fig.|3|it can be a crucial point how strongly 
a measurement of the line-of-sight velocity dispersion is contami- 
nated by unbound stars. The observational values can highly over- 
estimate the real central value. Unbound stars have a different ve- 
locity distribution than the bound ones. They are either travelling in 
front or behind the object (seen along the trajectory; this is shown in 
Fig-EJ and are faster or slower. But the velocity distribution of these 
stars which are still located in the vicinity of the object and stem 
from one central passage, peak on either side of the Gaussian distri- 
bution of the bound stars, leading to an effective broadening of the 
measured velocity distribution. Fitting just a single Gaussian will 
therefore lead to an enhanced measured velocity dispersion and an 
overestimation of the central velocity dispersion value of the satel- 
lite. This is demonstrated in Fig.|5|where the velocity distribution 
of all stars is shown together with the best fitting single Gaussian 
and the fitting curve of a triple Gaussian which takes the unbound 
stars into account. The measured values for the central line-of-sight 
velocity dispersion are shown in Table|5| 

Clearly this effect is strongest in the plane of the orbit and is 
not visible perpendicular to it. This can be seen in Fig. [5] where 
the symbol for the z— axis- value in all three panels shows the same 
value as the direct measurement. Actually for strongly disturbed 
systems (i.e. high mass-loss) the values in the z-direction are below 



Table 3. Central velocity dispersion along the x-axis derived by different 
methods. Shown are the values of the same three simulations plotted in 
Figs. |5]to|5]now labelled 'dissolved', 'enhanced' and 'massive', respec- 
tively. The rows show the different values of the central velocity dispersion. 
First row is the direct measurement of the central dispersion, 'single' de- 
notes the central dispersion derived from the fit of a single Gaussian to the 
data and 'triple (c)' denotes the value of a fit using three Gaussians where 
the central one (c) is used to compute the velocity dispersion of the bound 
object. 



[km/s] 


dissolved 


enhanced 


massive 


true actual 


10.6 ±3.1 


9.49 ± 0.05 


35.7 ±0.2 


single 


44.6 ±0.3 


15.90 ±0.30 


38.4 ±0.1 


triple (c) 


13.3 ±0.1 


12.35 ± 0.04 


35.3 ±0.3 



the line because we do not account for the change of the constant 
A which should increase for lower masses. But in any random ori- 
entation of the object with respect to the observer the effect of the 
enhanced M/M-ratio should at least be partly visible. 

Summing up our results we state the following: Satellites 
which are out of virial equilibrium, i.e. in the state of dissolution 
show enhanced virial masses and therefore the real mass content 
is overestimated. Figure [6] shows in the first panel the M/M-ratio 
plotted against the ratio of final to initial bound mass of the ob- 
ject. The dividing dashed line separates objects which have lost 
more than 50 per cent of their mass during the central passage and 
are already dissolved or are not likely to survive the next passage, 
from objects which are stable for several more passages through 
the centre. While on the left side the derived M/M-ratios (using 
the observational method) can climb up to very high values, the 
stable objects show only slight enhancements of the ratio if at all. 
But in the second panel one already sees that even if the object 
itself is in virial equilibrium (M v i r /Mb » 1.0) there are satel- 
lites which show a broadening of the velocity distribution leading 
to an enhanced mass-ratio if measured the observational way. The 
third panel finally shows an enlargement of this area and one finds 
M/M-ratios overestimated by up to a factor of five. Surviving ob- 
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Figure 6. Summary of the results derived using the observational method (Sect. l3~2*l . Af bs denotes the virial mass derived by the observational method, 
i.e. fitting a single Gaussian to the velocity distribution, Aft denotes the real bound mass of the object, M; n ; is the initial bound mass and Af v ; r is the virial 
mass derived by the analytical method (Sect. l3.il . Measurements along the x-axis are three-pointed stars (online: red), along the y-axis are crosses (online: 
green) and z-axis measurements are represented by five-pointed stars (online: blue). Simulations with initial parameters similar to the UCDs in Virgo are 
represented by filled symbols. First panel shows the derived mass-ratio as a function of the final/initial bound-mass-ratio. Satellites with final masses larger 
than 50 per cent of the initial mass (right of the dashed line) are deemed stable against immediate tidal destruction. Objects left of the dividing line are out 
of equilibrium and therefore have to show enhanced M/M-ratios. But there are also data points which are larger than one on the right of the dividing line. 
Second panel shows the same results but now the observational M/M-ratio is plotted against the analytical ratio. All data-points except the measurements 
perpendicular to the orbit are above the 'analytical' method (dashed line shows the line where the values of the observational method are equal to the direct 
method), z— values are below the dashed line because we do not correct the constant A for objects with high mass-loss. A clump of data-points is visible 
around M v - lr -/M^ = 1 which points to overestimated masses. The third panel finally is a magnification of this region. One clearly sees that the observational 
method can lead to overestimates of the true mass by a factor of up to five. Horizontal line shows a ratio of 1 with dashed lines denoting an error range of 
20 per cent. 



jects which show an enhancement in their M/M-ratio, if measured 
with this observational method, are already marked in Fig. |2| with 
a small 'e'. As one can see there is a range of critical distances 
(D > 100 pc and D < 1 kpc) to the centre of the host galaxy 
where an UCD like the ones found in Virgo could show an overes- 
timated mass-to-light ratio. 



3.3 Observability 

In the previous section we claimed that with the 'observational' 
method the mass-to-light ratios of UCDs could be overestimated, 
because the velocity distribution is not Gaussian any more but 'con- 
taminated' by unbound stars around the object. In this section we 
show that there is almost no chance for an observer to detect this 
'non-Gaussianity' of the velocity distribution. 

When measuring a spectrum of a distant object, observers 
have to deal with two major shortcomings. First there is the intrin- 
sic line-width produced mainly by the instrument. State-of-the-art 
instruments like UVES or Flames can reduce this line-width down 
to about 2 km s _1 . But the observations of the UCDs in Virgo were 
made with an instrumental line-width of about 25 kms" . 

The second effect an observer has to take into account is noise 
in the spectrum. 

Taking the velocity distribution from the enhanced M/M-ratio 
simulation we fold it with a Gaussian of the width <7i to mimic the 
instrumental line- width and determine at which line- width the 'fea- 
tures' of our distribution are washed out. This happens at an instru- 
mental line-width of <7i = 7.5 kms -1 (as shown in Fig.[7). Then 
we take the best state-of-the-art line- width of a\ — 2 kms -1 and 
add random white noise to the distribution until again the 'features' 
are almost invisible again. This happens already at a signal-to-noise 
ratio of 20 (see Fig. flower left panel). In the final panel of Fig.0 



we fold our distribution with the line-width of the observations of 
the Virgo-UCDs and add noise to mimic the same S/N-ratio as in 
the observations. There is no deviation from Gaussiantity visible 
any more. 



4 DISCUSSION & CONCLUSION 

We have shown with our models that dwarf satellites around a giant 
elliptical galaxy like M87 can have an enhanced mass-to-light ratio 
due to close passages to the centre of the host galaxy. While very 
close passages lead to the destruction of the satellite there are orbits 
which allow for enough 'damage' to the satellite to enhance the 
mass-ratio (measured the same way an observer would do) without 
completely destroying the object. 

The amount of destruction is larger if the satellite is less mas- 
sive and less concentrated. But the loss of about 20 per cent of the 
initial mass is enough to have the object surviving several more 
close passages and to mimic an enhanced mass-to-light ratio. For 
models comparable to the UCDs in Virgo and Fornax the range of 
possible minimum distances during central passages is about 100 to 
1000 pc. All passages closer to the centre lead to complete destruc- 
tion and all passages further away show no measurable effect at all, 
except for a mass-loss of the order of a few per cent. We there- 
fore conclude that the enhanced M/L-ratios measured for the Virgo 
UCDs by Hasega n et alj <2005l) may not be due to dark matter. 

A more general result of our study is the discrepancy between 
the derived virial masses if one has access to the correct proper- 
ties of the satellites compared to the virial masses derived the way 
an observer would measure. While very massive objects which are 
almost unaffected by tidal forces show the same results within the 
uncertainties one has to be careful if objects are less massive and 
are surrounded by a cloud of tidally stripped stars. These stars are 
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Figure 7. The top left panel shows the original velocity distribution. At the 
top-right panel this distribution was folded with a Gaussian of the width 
2 kms -1 which mimics the best intrinsic instrumental line- width available 
today. The middle panels show the same distribution but folded with a Gaus- 
sian of the width 7.5 kms - (left) and 25 kms -1 (right). The left panel 
shows the transition when the non-Gaussian features of the original distri- 
bution start to disappear. In the bottom panels we also added white noise to 
the folded distribution to mimic the noise an observer has in the spectrum. 
At the left panel we took the folded distribution with the best available line- 
width again and added noise to reach a S/N of 20. This is again the transi- 
tion where the features start to disappear. With more noise the features are 
not visible any more. Clearly with higher intrinsic line-widths the features 
disappear at lower noise-levels. The right panel mimics the properties of 
the spectra of the UCDs in Virgo i Haseaan et al. 2005). There the intrinsic 
line-width was 25 km s -1 and the S/N-ratio was 10. 



either faster or slower in the mean but their mean values are not too 
different to the bulk velocity of the bound stars to clearly disentan- 
gle the 'populations' (populations in the sense of within, in front or 
behind the object). Especially if the broadening of a spectral line 
is estimated by fitting the template line folded with a Gaussian for 
the instrumental line-width and a single Gaussian for the velocity 
distribution can the deduced central velocity dispersion be too high 
thus leading to a mass-to-light ratio that is too large by up to a factor 
of ten. 

Effects like this have to be given serious consideration when 



measuring velocity dispersions of faint and distant objects. 
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